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Einstein's equations of gravitation are not invariant under geodesic mappings, i. e. under a certain class of 
mappings of the Christoffel symbols and the metric tensor which leave the geodesic equations in a given 
coordinate system invariant. A theory in which geodesic mappings play the role of gauge transformations 
is considered. 
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1 Introduction 

If we start with Einstein's beautiful hypothesis that test particles in gravitational field move along geodesic 
lines of some Riemannian spaces V4, it is natural to expect that the differential equations for finding the 
metric tensor g a p{x) for a given distribution of matter also should be invariant under any transforma- 
tions at which the geodesic equations remain invariant. However, the geodesic equations are invariant not 
only under arbitrary transformation of coordinates (it is rather obvious) but also under geodesic mappings 
Pg 7 (x) — ► Tp^ix) of the Christoffel symbols in any fixed coordinate system [1]- [4]. 

If r^ 7 are Christoffel symbols in some coordinate system (x°, x 1 , x 2 , x 3 ), and we use coordinate time 
t = x a jc (c is speed of light) as a parameter along geodesic lines, then the differential equations of a 
geodesic line are of the form 

x a + (r^ 7 - c- 1 rg 7 i a )i /9 ±i' = 0, (l) 

where x a = dx/dt, 'x a = dx a /dt. It easily to verify that these equations are invariant under the mapping 

r^Or) = T a Pl {x) + 6% cf> 7 (x) + 5™<f> (x), (2) 

where <fi a (x) is a continuously differentiable vector field. In Riemannian space-time <fi a (x) can be ex- 
pressed in terms of determinants g and ~g of the metric tensors before and after the geodesic mapping of 
space-time V into V as follows: 



1 



F 7 

OL-y 



r' 



2(n + 1) dx a 



In 



(3) 



Such mappings of the Christoffel symbols in a given coordinate system induce some transformations 
not only of the curvature and Ricci tensors but also transformations g a p — * ~g af3 of the metric tensor which 
are obtained by solving of the partial differential equation 



9/3 T , a ( x ) = 2 <t><*{x)g M {x) + <j>p(x)g ia {x) + 4> 7 (x)g a p(x), 



(4) 
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where the semicolon denotes a co variant derivative with respect to the metric in V4. 

It is clearly that all solutions of these equation are, in principle, equivalent physically. However Ein- 
stein's equations even in vacuum are not invariant with respect to geodesic mappings [5]. 

It follows from that the components Tq = r oo . Therefore, in Newtonian limit geodesic-invariance 
is not an essential fact. However this cannot be said about the relativistic case. For this reason is of interest 
to explore a theory in which the gravitation equations, as well as the equations of motion of test particles, 
are geodesic- invariant i.e. in which geodesic invariance plays the role of gauge invariance. 

In Sect. 2 a geodesic-invariant generalization of Einstein's equations is considered. It turns out that the 
simplest equations of such a kind are some bimetric equations. In Sects 3 to 7, proceeding from new sight 
at Poincare old fundamental ideas, we argue that the both spaces, Riemannian (with the curvature other 
than zero) and at, can have physical sense depending on which kind of reference frame is used. In Sects. 
8, 9 we find the spherically-symmetric solution of the above equations. Finally, in Sect. 10 it is shown 
that features of the gravitational force resulting from this solution, well explain properties of the Hubble 
diagram obtained from the latest observations. 



2 Generalization of Einstein's equations 

The simplest way to find a gauge-invariant generalization of Einstein's equations is to consider compo- 
nents T^g of the affine connection in n — 4 -dimensional manifold M.4 as 4-components of some affine 
connection Tg C of n + 1 dimensionall manifold in which transformations (0 are consequences of 5th 
coordinate transformations [2]: 



-4 4 
X = X 



foWdx'P , (5) 







where <t>p(x') is a gradient vector field. These transformations together with 4-transformations 

x a = x a (x°,x 1 ,x 2 ,x 3 ) (6) 

form some acceptable coordinate transformations in A4s . In one can define an affine connection by 
the functions which are transformed under (0 and (0 as follows: 

-A _ ( F dx E dx D d 2 X F \ CfX A 

The equations for 4-components: 

=« _ ( vF dx E dx B d 2 x F \ dx<* 



under transformations (0 takes the form (0 if Christoffel's symbols satisfy the conditions : r? 4 = T" 7 = 
6° and r Q 4 = 0. These conditions are invariant under coordinate transformations (O and ©. The rest 
components of Tg C should satisfy the conditions which also are obtained from consideration of their 
transformation properties under (01. The components Tq 7 = in (O do not depend on T FD . This allows 

to set Tq 7 = in any coordinate system. The component T 44 are equal to T| 4 , which allows to set 
T 44 = 1 in any coordinate system. The components T* « are transformed as a tensor under coordinate 
transformations (O, and as 



—4 



r a/3 = r a/3 - i>a/3 (9) 

under transformations (0 of the 5th coordinate where ip a p = ij) a -p — ipaipp, and tpa-,/3 is a covariant 
derivative in V4. 
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The object 

T^bcd = dTi c /dx D dTi D /dx c + T^T^c - Ti w T A J D (10) 

is a curvature tensor in M.§ . In our case the components lZi a = 7l a 4 an d the component 7^.44 of the 
contracted tensor TZ bc = 'Rbac 816 e( J ua l to zero identically. For this reason TZp 7 is the 4- tensor 

ft/3 7 = Rm + (n - 1)1^, (11) 

where i?^ 7 is the Ricci tensor. 

The components of the Ricci tensor are transformed under (O as follows: 

Rap = RaP — lpa/3- (12) 

On account of (O the object IZp^ is gauge (i.e geodesic) invariant. For this reason, the simplest geodesic- 
invariant generalization of vacuum Einstein's equations is of the form 

R 01 + (n - l)r£ 7 = 0, (13) 

where r^ 7 is a tensor field which is transformed under geodesic mappings according to (0. 

Now the problem is to find some relevant tensor P| . Such object must satisfy the following require- 
ments: 

1 . It should transform under (0 according to (0, 

2. It should transform under © as a 4-tensor, 

3. It should lead to such generalization of Einstein's equations which solutions can noticeably differ 
from the solutions of the last equations only in strong gravitation field as in weak and moderately strong 
field solutions of Einstein equations are in good agreement with observations. 

The first requirement means that should be formed from geometrical objects of space-time V4, 
and cannot be an additional field. Because the equations of motion of free particles contains connection 
coefficients (and not metric) it should be formed from components of H . 

The simplest object T^ 7 which can be formed by components of Christoffel symbols and which has the 
requirement transformations properties under our gauge transformations is 

T ^=^ T ^- r ^-jdw TaT0 ' (14) 

where T a = . In this case IZp^ is an object which is formed by gauge-invariant Thomas symbols [2] 



(15) 



the same way as the Ricci tensor is formed out of the Christoffel symbols. As II" = 0, equations (TT3T > 
take the form 

v 7 n« - n^nf 7 = o (16) 

where V a denotes a covariant derivative in Minkowski space-time. 

However, the problem is that TZp^ is not a tensor with respect to general coordinate transformations 
because r Q (x) is not a vector field. 

This a new problem evidently cannot be solved within the framework of General Relativity because by 
using components T^ 7 it is impossible to form a tensor of the second rank in V4 which is invariant under 
geodesic transformations. 

However a satisfactory solution of this problem can be found if the Riemannian space-time V4 coexists 
in the manifold M4 with Minkowski space-time E4 which also has a physical meaning. From the point of 
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view of such bimetric structure of space-time, in Minkowski space-time g a p(x) is simpy some tensor field. 
The Christoffel symbols T^p can be considered as components of the tensor = Pg — Tg^ in Cartesian 
coordinates, where Ta is the Christoffel symbols in Minkowski space-time. Tensor Dg is formed by 



substitution of ordinary derivatives in by the covariant ones0because both the definitions of the tensor 

o 

coincide in Cartesian coordinates. The object T a is components of the vector D a = = — T* 3 a p 
in Cartesian coordinates. 

In such bimetric manifold we can set 



a0 



}a,f3 — T^c, Qj — Q a Qfj — Q a ;p — Q a Qp, 



where 



n+1 



-D r 







2(n+ 1) dx a 



In 



The definition (fl4l i is the tensor (T7\ in Cartesian coordinates. 

Similarly to it, instead of Thomas's symbols (TT~5T > we have the gauge-invariant tensor 



b 1 - n 1 

D aP 



It can be written also as 



51 Qp + 6} Q a 



DO; Y\a 



n 



Pi 



where 



o7 



o7 



n 



a/3 



r a/3 -(n + i)- x 



are Thomas symbols for £4 in the used coordinate system. 

Therefore, Eqs. ( fTol l can be considered as a particular case of the more general equations 



V Q Bg 7 - B e gs B t 



0. 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



where V Q denotes a covariant derivative in E4,. 

The above can be generalized to the case of matter presence. It is naturally to consider the tensor g a p 
as 4-components of some metric tensor gAB in the manifold M.5 : 



gAB 



gap g a i 

<74a <744 



The components g a p are transformed under (O as follows: 

9 a p = Sap + ga4> a (t)p + gia4>p + gip4>a 

g ai = 9ai + g<LA<f>a 
9a4 = .944- 

Transformations of the components T a under geodesic mappings are given by 

T a = T a + (n + l)<t> a , 



(23) 



(24) 
(25) 
(26) 



(27) 



There is an analogy here with the formalism of the bimetric Rosen theory [33]. 
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which coincides with the transformations of the components g a ^ if 544 
assume that0 

9AB= { Qa n+1 )■ 
Then there is a gauge-invariant tensor 

g a p = g a /3 -{n+ l)~ 1 Q a Qp- (29) 
Thus, finally, the gauge-invariant generalization of Einstein's equations with matter are of the form 

V Q S^ 7 - Bp S B^ = k ^T a p - -g a pT^j , (30) 

where k = Sttj/c 4 , 7 is the gravitational constant, c is the speed of light, T a p is the matter energy- 
momentum tensor. At the gauge conditions Q a = equations ( f30b coincides with Einstein's equations. 

Of course, these bimetric equations may be true if and only if both the space-times, V4 and £4, has 
some physical meaning. How these two physical space-time can coexist? Attempts to use flat space-time 
simultaneously with curved in general relativity were undertaken repeatedly. Here it should be mentioned 
the well known Rosen's bimetric theory [6] which leads to non-Einstein equations for g a p (x) , and inter- 
pretations of the Einstein's equations as the equations for g a p considered as some function of a tensor field 
ipa/3 which describes gravitational field in Minkowski space-time [7], [8], [9], [10]. The need to go beyond 
Riemannian manifold has noted been also in [11]. 

It should be noted that at any concrete expressions of g a p in terms of tp a p, which are used in some 
papers, the Einstein equations, considered as a function of ip a p, cannot be invariant with respect to pos- 
sible gauge transformations of the field ipafi- In the paper [12] it was shown that the requirements of the 
invariance of the any equations for g a p with respect to the possible gauge transformations of ip a p leads 
uniquely to geodesic invariance of gravitational equation, and the simplest generalization of Einstein's 
vacuum equations are equations ( fT3l as obtained above. 

However bimetric nature of physical space-time belongs probably to one of the deepest problems of 
physics. In Sec. 3 to 7 we show how both the kind of space-time can coexist in the theory. 



1. For this reason we can 



(28) 



3 Relativity of space-time 

At the beginning of the 20th century H. Poincare convincingly showed that it makes no sense to assert that 
one or other geometry of physical space is true. Only an aggregate " geometry + measuring instruments" 
has a physical, meaning verifiable by experiment. Einstein recognized that Poincare reasoning is justified. 
His theory of gravity does not reject Poincare's ideas. It only demonstrates relativity of space and time 
geometry with respect to distribution of matter. Poincare's ideas indicate to another property of physical 
reality — relativity of the geometry in relation to measuring instruments, dependency on their properties. 
These ideas, apparently, have never been realized in physics and, evidently, remain until now a subject for 
discussion of philosophers. The success of General Relativity almost convinced us that physical space- 
time in the presence of matter is Riemannian and that this fact does not depend on properties of measuring 
instruments. Attempts to describe gravity in flat space-time have not obtained recognition and Poincare's 
ideas have proved to be almost forgotten for physics. 

However, are Poincare's ideas really only a thing of such little use as conventionalism, as it is usually 
believed? In fact, a choice of certain properties of measuring instruments is nothing more than the choice 
of some frame of reference, which is just such a physical device by means of which we test properties of 
space-time. (Of course, we must understand the distinction between a frame of reference as a physical de- 
vice, and a coordinate system as only a means of parameterization of events in space-time.) For this reason 



It means that there is some relationship between Tg C and g^B distinct from such relationship in Riemannian geometry. 
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Poincare's reasoning about interdependency between properties of space-time and measuring instruments 
should be understood as the existence of a connection between geometry of space-time and properties of 
the employed reference frame . 

The existence of such a connection can be shown by an analysis of a simple and well-known example 
considering it from a new point of view. Disregarding the rotation of the Earth, a reference frame, rigidly 
connected with the Earth surface, can be considered as an inertial frame (IFR). An observer, located in 
this frame, can describe the motion of freely falling identical point masses as taking place in Minkowski 
space-time under the action of a certain force field. However, consider another observer, who is located 
in a frame of reference, the reference body of which is formed by these freely falling particles. Such a 
reference frame can be named the proper frame of reference (PFR) for the given force field. Let us assume 
that the observer is deprived of the possibility of seeing the Earth and stars. This observer does not feel 
the presence of the force field in any place of his frame. Therefore, if he proceeds from the relativity of 
space-time in the Berkeley-Leibnitz-Mach-Poincare (BLMP) sense [13] - [16], then from his viewpoint 
accelerations of the point masses, forming the reference body of his frame, in his physical space must be 
equal to zero. However, instead of this, he observes a change in distances between these point masses in 
time. How can he explain this fact? Evidently, the only reasonable explanation for him is the interpretation 
of this observed phenomenon as a manifestation of the deviation of geodesic lines in some Riemannian 
space-time of a nonzero curvature. Thus, if the first observer, located in the IFR, can postulate that space- 
time is flat, the second observer, located in a PFR of the force field, who proceeds from relativity of space 
and time in the BLMP sense, already in the Newtonian approximation is forced to consider space-time as 
Riemannian with curvature other than zero. 

In next Section it is shown that a similar dual description is possible also for the motion of a perfect 
isentropic fluid. 



4 Space-time metric form in NIFRs 

It is possible to arrive at some quantitative relationship between properties of the used frames of reference 
and the local metric properties of space-time in PRFs examining the problem from more general positions. 

Presently we do not know how the space-time geometry in an inertial frame of reference (IFR) is related 
with the frame properties. Under the circumstances it can simply be postulated, according to Special 
Relativity, that space-time in IFRs is pseudo-Euclidean. Proceeding from it, one can find the line element 
of space-time in non-inertial frames of reference (NIFR) from the viewpoint of observers located in the 
NIFRs and proceeded from relativity of space and time in the BLMP meaning. 

By a non-inertial frame of reference we mean the frame, the body of reference of which is formed by 
point masses moving in an IFR under the effect of a given force field. 

The reference body (RB) of a reference frame is supposed to be formed by identical point masses m. If 
an observer in the reference frame is at rest , his world line coincides with the world line of some point of 
the reference body. It is obvious for such an observer in an IFR that the accelerations of the point masses 
forming his reference body are equal to zero. This fact takes place in relativistic meaning, too. That is, if 
the line element of space-time in the IFR is denoted by da and u a = dx a /da is the field of 4- velocity of 
the point masses forming the reference body, then the absolute derivative of u a is equal to zero: 

Du a /da = 0. (31) 

(We mean that an arbitrary coordinate system is used.) 

Does it occur for an observer in NIFRs ? That is, if the differential metric form of space-time in a NIFR 
is denoted by ds, does the 4- velocity vector ( a = dx a /ds of the point-masses forming the reference body 
of the NIFR satisfy the equation 

DC/ds = ? (32) 



We use notations and definitions, following the Landau and Lifshitz book [22]. 
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If the space-time is absolute, equation (l32l holds for only ds = da. However, if space and time are relative 
in the BLMP sense, then for both observers, located in some IFR and NIFR, the motion of the point masses, 
forming their reference bodies, which are kinematically equivalent, must be dynamically equivalent, too 
(both in non-relativistic and relativistic sense). Any observer in the NIFR, isolated from the external world 
and proceeded from relativity of space-time in BLMP meaning, consider points of the reference body as 
the ones of his physical space, and space of events as his space-time. Therefore, from his viewpoint point 
masses forming the reference body of his frame are not under action of any forces (the same as for the 
observer in IFR), and their 4-velocity must be equal to zero. In other words, since for the observer in the 
IFR world lines of the reference body are, according to PTT ). some geodesic lines, for the observer in the 
NIFR the world lines of the of his BR also must be geodesic lines in his space-time, which can be expressed 
by (QU). 

The equation d32l uniquely determines the fundamental metric form in NIFRs. Indeed, the differential 
equations of these world lines are at the same time Lagrange equations describing in Minkowski space-time 
the motion of the point masses forming the reference bodies of the NIFR. The last equations can be obtained 
from a Lagrange action S by the principle of the least action. Therefore, the equations of the geodesic lines 
can be obtained from a differential metric form ds = k dS, where k is constant, dS = C(x,x)dt, and 
C(x, x) is a Lagrange function describing in Minkowski space-time the motion of identical point masses 
m forming the body reference of the NIFR. . The constant k is equal to — (mc) _1 , as it follows from the 
analysis of the case when the frame of reference is inertial. 

Thus, if we proceed from relativity of space and time in the BLMP sense, then the line element of 
space-time in NIFRs can be expected to have the following form 

ds = — (mc)^ 1 dS(x,dx). (33) 

Therefore, properties of space-time in NIFRs are entirely determined by properties of used frames in 
accordance with the BLMP idea of relativity of space and time, which has been noted for the first time 
in [19]. 

Consider two examples of NIFRs. 

1 . The reference body is formed by noninteracting electric charges, moving in a constant homogeneous 
electric field £ . The motion of the charges in an IFR is described in the Cartesian coordinates system by a 
Lagrangian [22] 

L = -mc 2 (1 - v 2 /c 2 ) 1/2 + £ ex, (34) 

where v is the speed of a particle. According to d33l the space-time metric differential form in this frame 
is given by 

ds = da — (wx/c 2 )dx°, (35) 

where 

da = [r, a pdx a dx P ] 1/2 

is the differential metric form of the Minkowski space-time in the IFR in the coordinate system being used, 
and w = e£/mis the acceleration of the charges. 

2. The reference body consists of noninteracting electric charges in a constant homogeneous magnetic 
field H directed along the axis z. The Lagrangian describing the motion of the particles can be written as 
follows [22]: 

L = -mc 2 (l - v 2 /c 2 ) 1/2 - {mn /2)(xy - xy), (36) 



where x = dx/dt, y = dx/dt, and SIq = eH/2mc. 
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The points of such a system rotate in the plane xy around the axis z with the angular frequency 

W = ft [l + (tt r/c) 2 r 1/2 , (37) 

where r = (x 2 + y 2 ) 1 ^ 2 - The linear velocity of the BR points tends to c when r — > oo. 
For the given NIFR 

ds = da+ /2c) (ydx - xdy). (38) 
In the above NIFRs ds is of the form 

ds = J-{x, dx), (39) 

where T{x, dx) = da + f a (x)dx a , f a is a vector field. Therefore, !F(x, dx) is a homogeneous function 
of the first degree in dx a , that is !F(x,ky) = hF{x : y). Thus, the space-time in the above NIFRs is 
Finslerian [23]. 

By using the identities 

T(x y) = 9F{x ' y) r- g!^M) =Q ( 40) 

the function !F 2 (x, y) can be written as 

T{x,yf = G aP {x,y)y a y l3 1 (41) 

where 

g ^tf) = 2 (42) 

is an analog of the fundamentall metric tensor of Riemannian space-time. Consequently, the modulus of 
the vector y a in the point a;" is \\y\\ =J 7 (x,y) = [G a p(x,y)y a yP] 1 / 2 . 
For ds of the form d39b we have 

G a p{x, y) = rj af 3 + f a fp + {r]af3 ~ T a Tf3)f a T a + f a T + fpr a , (43) 

where r a — dx a /da. 

A covariant vector in the Finslerian space-time, therefore, can be defined as 

L = T{x, V)^^ = G a p(x, y) y a . (44) 

The orthogonality of vectors y a and y" can be defined by the equality y a yf = 0. It must be noted that 
the orthogonality of two vectors is not symmetric. 

3. The reference body is formed by particles of an isentropic flow of the perfect fluid. 

It is well known that besides of traditional continual description, a perfect fluid can be considered 
as a collection of a finite number of identical macroscopic small particles which are under influence of 
interparticles forces which mimic the effect of pressure, viscosity etc [17] , [18]. In particular, the fluid 
velocity in a given point is simply velocity of the particle being in this point. At such description the 
motion of the fluid particles is governed by solutions of ordinary differential equations of Newtonian or 
relativistic dynamics. 

The motion of the fluid particles in an IFR is described by Lagrangian [24] 



dx a dx 



1/2 



- mc l G ^^r) dA ' (45) 
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where A is a parameter along 4-path of particles, w is enthalpy per unit volume, G a p = x 2 rj a p , x = 
w/pc 2 , p = ran, m is the mass of the particles, n is the particles number density, and r) a p is the metric 
tensor in Minkowski space-time. According to d33l l the line element of space-time in this frame is given by 

ds 2 = G aP dx a dx (3 . (46) 

The covariant derivative of tensor G a p is equal to zero. Therefore, space-time in such NIFR is Riemannian 
with the curvature other than zero. The last description is equivalent to the one by the Lagrangian d45l l [24]. 

Thus, an observer in some laboratory (inertial) frame of reference can describe a fluid motion as what is 
going on in Minkowski space-time by a Lagrangian. However, the observer which is located in a comoving 
frame of reference cannot feel presence of the force field arising from a pressure gradient since he moves 
along a geodesic line of some Riemannian space-time. Therefore, if he is isolated from external world, 
he is forced to explain relative motions of nearby particles of the fluid as a manifestation of deviation of 
geodesic lines due to space-time curvature in this frame of reference. 



5 3+1 decomposition of space-time in NIFRs and Sagnac effect 

For the 3 + 1 decomposition of space-time in non-inertial frames of reference into physical 3-space and 
time we use a covariant method which is a Finslerian generalization of the well-known method in General 
Relativity [25], [26]. 

An ideal clock is a local periodic process measuring the length of its own world line on a certain scale. 
For an observer in a NIFR the direction of physical time in a point x a is given by the vector of the 4-velocity 
( a = dx a jds of the given point of the NIFR reference body. 

By analogy with the case of Riemannian geometry, an arbitrary 4-vector £ a in the point x a can be 
represented as follows: 

c a =r (47) 

— OL 

In this equation p is a function of x a , and £ is a 3-spatial component of the vector It is supposed to 
be satisfied to Finslerian's orthogonality condition [23] to the time direction which we define as 

C«r = 0, (48) 
where ( a is the covariant components of the vector £ a : 

Since T{x, £) = 1, this vector is 

c ;= Va ^ al/2 +f a . (so) 



By multiplying d47b with the vector £ Q we find that (3 = C Q £ Q an d 

'/3 



where Hp = C, a Qp — Sg is an operator of spatial projection, and 6% is the Kronecker delta. 
For the vector £ Q = dx a d47T > yields 

dx a = dx~ a + cdT f ( a , (52) 
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where dx is the spatial components of the vector dx a , and 

dT f = c -1 C a dx Q (53) 

is the time element between events in the points x a and x a + dx a in the NIFR. 

Metric form d39t and the spatial projection of the vector dx a leads to the following covariant form of 
the spatial element in NIFRs 

dLf = (-r^d^dz 3 ) 172 + /a dx~ a . (54) 
This covariant equation is the simplest and clearest in the "comoving" coordinate systems, in which 

C=e5%. (55) 

It follows from the equality J-(x, Q = 1 that e = (i]Iq 2 + fa) 1 - 

In this coordinate system for an observer at rest the time element dTf = (?7qq 2 + fo)dx°, dx 1 = dx 1 , 
and dx° = 0. Therefore, 

dLf = {~-T lik dx i dx k ) 1/2 + f l dx i = dL(l + fck 1 ), (56) 

where dL = (—iiikdx l dx k ) 1 / 2 is the Euclidean spatial element and k % = dx 1 jdl is the unit direction 
vector. 

Of course, the above definitions of dTf and dL f have physical meaning only in the space-time region 
where they are real and positive quantities. In this respect, space-time in NIFRs has more complicated 
structure as compared with IFRs, and depends on properties of the considered NIFR. For example, space- 
time in the first NIFR in Sec. 4 is time-like only if the interval da is time-like and the distance \x\ of 
charges from the coordinate origin satisfy the condition \x\ < c 2 /w. 

The Finslerian geometry of space-time in NIFRs allows us to give a natural explanation of the well 
known Sagnac effect, i.e. of a phase shift in the interference of two oppositely directed coherent light 
beams on a rotating disk [27]. It is usually, for relativistic explanation of this effect the motion of light in a 
NIFR is considered as a relative motion in the Pseudo-Euclidean space-time of some IFR [28]. However, 
for an observer located in a rotating frame and isolated from the IFR (in a "black box") , who proceeds from 
the notion of space and time relativity in BLMP sense , the observed anisotropy in time of light propagation 
clockwise and counterwise is not a trivial effect. It must have some "internal" physical explanation. 

Consider a rigid disk rotating in the plane xy with a constant angular velocity il around the axis z. Let 
r and 8 be the coordinates, defined by the equations 

x = r cos ip, y = r sin ip, tp = 9 + fit. (57) 
In the coordinate system (r, 6, z, t) the differential metrical form ds is of the form 

ds = da- {VLr 2 /2c) dO - (n 2 r 2 /2c 2 )dx° , (58) 
where da is the pseudo - Euclidean metric form: 

da 2 = (1 - n 2 r 2 /c 2 ) dx° 2 - dr 2 - r 2 d6 2 - 2(r 2 n/c)d6 dx° - dz 2 . (59) 

At HR/c <C 1 this reference frame is identical to the NIFR described in example 2 of Sec. 4, and the 
used coordinate system on the disk is "comoving", where d55b holds. 

It follows from d56b that the spatial element in the NIFR is anisotropic. We will show that the speed of 
light in the noninertial frame of reference is anisotropic, too. Along photon paths the equality r\ a ^dx a dx^ = 
holds. By using d52l , this equation can be written in terms of notions of the NIFR as 

ri^dx^dx 13 + 2cdT f Tj a pCdx P + c 2 dTjr) af3 C( = 0, (60) 
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where the 4-vector dx is orthogonal to ( a , i.e. 



fa 



dx =0. (61) 



For this reason the preceding equation takes the form 

W^ h v% ~2cf a v° h (r la0 CC ) 1/2 + c 2 r, afS Ct = 0, (62) 

where v™ h = dx" /dTf is a photon velocity as measured by an observer in the NIFR. 

The covariant equation d62l become more clearly in the coordinate system where the conditions ( f55l > 
are satisfied. 

Indeed, it follows from the equality F(x,() = 1 and from 051) and (T58l that the quantity (rjapC^C 13 ) 1 ^ 2 = 
1 — f a ( a = 1 — /2C 2 — /oC° = 1 w i m accuracy up to Qr/c. Therefore, photon's velocity v p h in the NIFR 
satisfy the equation 

v 2 h + 2cv ph f i k i -c 2 = 0, (63) 

where 

2 dx 1 dx k 

v Ph = -mk WfWf . (64) 

Thus, the velocity of an outgoing photon is given up to accuracy Slr/c by the equation 

v ph = c(l-f i k i ). (65) 
The time of the photon motion around of the disk of the radius R is an integral of 

dTf = d ^l = ^{l + 2f i k% (66) 

Vph c 

which yields for clockwise and counterclockwise directions dTp+ = (dL/c)(l + 2/g) and dTp- = 
(dL/c)(l — 2/e), respectively. The difference in the time interval between light propagation on the rotating 
disk in two directions is AttR 2 Q,/c 2 , which gives the Signac phase shift [27]. Thus, the Signac effect for 
an isolated observer in the rotating frame can be treated as caused by the Finslerian metric of space-time 
in non-inertial frames of reference, see also [29]. 

6 Inertial Forces 

Let us show that the existence of inertial forces in NIFRs can be interpreted as a manifestation of a Finsle- 
rian connection of space-time in these frames. 

According to our main assumption in Sec. 2, the differential equations of the motion of the point 
masses in an IFR , forming the reference body of a NIFR, are equations of geodesic lines of space-time in 
this NIFR. These equations can be found from the variational principle 5 J ds = and are of the form 

d( a /ds + G a (x,() = 0, (67) 
where ( a is the 4-velocity of the point mass, the world line of which is x a = x a (s), and 

G a (x,Q =T a M CC +B$C P + C Zd{Z- 1 )/ds. (68) 
In these equations Pg are the Christoffel symbols in Minkowski space-time, 

Z- 1 = da/ds = (vapC( f3 ) 1/2 , B% = V aS Bsp 
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and 

B S0 = d.f, 3 /dx & - dfs/dx 13 . 

In Finslerian space-time a number of connections compatible with d6Tt can be defined [23]. In particu- 
lar, this equation can be interpreted in the sense that in NIFRs space-time the absolute derivative of a vector 
field £, a (x) along the world line x a = x a (s) is given by 

DC/ds = dC/ds + G%{x, dx/ds)f, (69) 

where 

G%(x,dx/ds) = B% + T^ ( dx" / /ds + ZdZ^/ds. 

Equations J69b define a connection of the Laugwitz type [23] in the space-time of a NIFR, which is 
nonlinear with respect to dx a . The change in the vector due to an infinitesimal parallel transport is 

d£ a = -G%{x,dx)£ p , (70) 

Consider a free motion of a particle of mass m in a NIFR. Since ds = da — f a dx a , equations of the 
motion are 

Du a /ds = B%u !i . (71) 

The equations of the motion ( f67b on a non-relativistic disk rotating in the plane xy about the axis z with 
an angular velocity £1 are take the form 

dC/dt + Q x f = 0, , (72) 

where r = {x, y, z} and the coordinates origin coincides with the disk center. The absolute derivative d69b 
of a vector £ is given by 

D£/dt = dg/dt-tix £ (73) 
and equations of the motion ( fTTl ) of the considered particle in the NIFR are 

Dv/dt = -Qxv, (74) 

where v = {xy, z}. 

Next, for the 4-velocity u a we have 

u a = u a +( a , (75) 

where u a is the spatial velocity of the particle in the NIFR. In the non-relativistic limit ( fTBI l can be written 
in the form 

v = ^+(, (76) 

where v is the "relative velocity" of the particle and £ is the velocities field of the points of the disk in the 
laboratory frame. Substituting ( 1761 ) for d74l . we find that 

Dlj/dt = -D(/dt-<rixv-nx(. (77) 

The value Dv / dt is an acceleration of the considered particle as measured by an observer in the NIFR. 
The velocities field ( of the disk points is given by 

( = O x f. (78) 
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Hence, along the particle path we have dQJdt = X v and 

D(/dt = d(/dt - (1 X ( = X (79) 
Therefore, finally, we find from d74b 

mDlJ/dt = -2m(Sl x H) - mQ X (f2 X r). (80) 

We have arrived at non-relativistic equations of motion of a particle in rotating frames [30]. The right-hand 
side of ( f80b is the expression for the Coriolis and centrifugal forces. 

Thus, in non-relativistic limit the Finslerian space-time in NIFRs manifests itself in the structure of the 
derivative of vectors with respect to time t , see also [29]. 

It is interesting to note that d73l l is considered sometimes in classical dynamics nominally [30] just to 
obtain the expression for inertial forces in NIFRs. 



7 Relativity of Inertia 

A clock, which is in a NIFR at rest, is unaffected by acceleration in space-time of this frame. A difference 
in the rate of an ideal clock in IFRs and NIFRs is a real consequence of a difference between the space-time 
metric in the IFR and NIFR. It is given by the factor Z = ds/ da. For a rotating disk of the radius R the 
factor Z = 1 — uj 2 R 2 /2c 2 which gives rise to the observed redshift in the well-known Pound-Rebka-Snider 
experiments. 

Consider another experimentally verifiable consequence of the above theory. Let p a = mc dx a /da be 
4-momentum of a particle in an IFR. From the point of view of an observer in a NIFR 

p a = p a +c- 1 E( a . (81) 

In this equations the spatial projection p a should be identified with the momentum and the quantity E with 
the energy of the particle. 

It is obvious that E = ( a p a - Therefore, the energy of the particle in the NIFR is 

E = mZc 2 l a u a , (82) 
where Z = ds/da = F(x, dx/da). For the particle at rest in the NIFR u a = ( a and we obtain 

E = mZc 2 . (83) 
Thus, the inertial mass mi of the particle in the NIFR is given by 

mi = Zm. (84) 

The quantity m coincides with the proportionality factor between the momentum and the velocity of a 
nonrelativistic particle in the NIFR. 

Since Z is the function of x a , the inertial mass m n in the NIFR is not a constant. For example, on a 
rotating disk we have 

m n = m/{l-n 2 r 2 /2c 2 ), (85) 

where SI is the angular velocity and r is the distance of the body from the disk centre. The difference 
between the inertial mass m eq of a body on the Earth equator and the mass mP° l of the same body on the 
pole is given by 

{m eq ~m pol )/m po1 = 1.2 • 10~ 12 . (86) 

Copyright line will be provided by the publisher 



16 



L. Verozut|: Geodesic-invariant equations of gravitation 



The dependence of the inertial mass of particles on the Earth longitude can be observed by the Mossbauer 
effect. Indeed, the change AA in a wave length A at the Compton scattering on particles of the masses m is 
proportional to m. If this value is measured for gamma-quanta with the help of the Mossbauer effect at a 
fixed scattering angle, then after transporting the measuring device from the longitude ip\ to the longitude 
(fi2 we obtain 

A \— 1 _ A \— 1 

Vl _ 1 V2 = 9 [cos 2 Vl - cos 2 ip 2 ] , (87) 

AA^j 

where 6 = 1.2- 1(T 12 . 



8 Gravitation in Inertial and Proper Reference Frames 

To implement the the bimetric description of gravity based on the relativity of space-time to the employed 
reference frame , suppose [7] that in Minkowski space-time gravitation can be described as a tensor field 
4>ap{x) of spin 2, for which the Lagrangian, describing the motion of a test particle of mass m is of the 
form 

L = -mc[g a p{^) x a x 13 ] 1 ' 2 , (88) 

where x a = dx a /dt and g a p is the symmetric tensor whose components are functions of tjjafS- 

Consider a frame of reference, the reference body of which is formed by identical point masses m 

moving in an IFR under the effect of the field ip a /3 (x). It is a proper frame of reference of the given field. 

An observer, located in the PFR at rest, moves along a geodesic line of his space-time. 
According to d33l the square of the line element in a PFR of the field ip a p is given by 

ds 2 = g a p{i>) dx a dx p . (89) 

Thus, according to what has been outlined in Sec. 4, space-time in PFRs is Riemannian V4 with the 
curvature other than zero. Viewed by an observer located in the IFR, the motion of the particles, forming 
the reference body of the PFR, is affected by the force field i/Ja/3- Let x l (t, x) be a set of paths of the 
motion of the depending on the parameter \- Then, for the observer located in the IFR the relative motion 
of a pair of particles from the set is described in non-relativistic limit by the differential equations [32] 

dt 2 dx l dx k ' ^ 

where n k = dx k jd\ and U is the gravitational potential. 

However, the observer in a PFR of this field will not feel the existence of the field since he moves in 
the space-time of the PFR along a geodesic line. The presence of the field ip a p will be displayed for him 
differently — as space-time curvature which manifests itself as a deviation of the world lines of nearby 
points of the reference body. 

For a quantitative description of this fact it is natural for him to use the Riemannian normal coordi- 
nates [3 In these coordinates spatial components of the deviation equations of geodesic lines are 

H+iWi fc =o, (9i) 

where i?Q fe0 are the components of the Riemann tensor. If the metric tensor does not depend on time t, then 

R i _ ri 9 goo 

Ron ~-9 g^- (92) 



^ This and the above consideration does not depend on the used coordinate system, it can be performed by a covariant method. 
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In the Newtonian limit these equations coincide with d90l ). 

Thus, in two frames of reference being used we have two different descriptions of particles motion — 
as moving under the action of a force field in Mankowski space-time, and as moving along the geodesic 
line in a Riemann space-time with the curvature other than zero. 

Of course, ( T33l > refers to any classical field. For instance, space-time in PFRs of an electromagnetic field 
is Finslerian. However, since ds, in this case, depends on the mass and charge of the particles forming the 
reference body, this fact is not of great significance. 

If we start from the Lagrangian (l88l for the motion of test particles, Einstein's equations cannot be con- 
sidered as equations for finding functions g a p (x) . The reason is that the gravitational equations for g a p (if)) 
should be some differential equations which are invariant under a certain group of gauge transformations 
ipap(x) — ► ipap (x) which are a consequence of the existence of "extra" components of the tensor ip a p. 
These transformations induce some transformations g a p{x) — > ~g a p{x). Equations of motion of particles 
resulting from the Lagrangian (l88l certainly should be invariant under such transformations of the tensor 
9af){%). But the equations of motion are, at the same times, differential equations of geodesic lines in 
Riemannian's space-time of a proper frame of reference. Consequently, if equations of gravitation do not 
contain functions ip a p explicitly, i.e. are differential equations for g a p(x), they should be invariant under 
such transformations of metric tensor of space-time in PRFs, and it takes place in any coordinate system at 
that. 



9 Spherically-Symmetric Gravitational Field 

Let us find a spherically-symmetric vacuum solution g a p(x) of bimetric equations ( f30b in Minkowski 
space-time E4. Such a solution, considered as a tensor field in E4, describes, in particular, the gravity of a 
point mass M for a remote observer in an IFR for which space-time is supposed to be pseudo-Euclidean. In 
this case, if the Lagrangian of a particle is invariant under the mapping t — » —t, in the spherical coordinates 
of the Minkowski space-time, it is of the form 

L = -mc[Ar 2 + B(6 2 + sin 2 9 ip 2 ) - c 2 C] 1/2 , (93) 

where A, B and C are the functions of the radial coordinate r. 
The associated line element of space-time in PFRs is 

ds 2 = Adr 2 + B[d9 2 + sin 2 6 <p 2 } - C dx° 2 \ . (94) 

Some additional conditions can be imposed on the tensor g a p because of the gauge (geodesic) invari- 
ance. In particular, under the conditions Q a = Eqs. d22T > are reduced to Einstein's vacuum equations 
Rap = 0.. Thus, the functions A, B and C can be found as a solution of the system of the differential 
equations : 

R a p = (95) 

and 

Q a = 0, (96) 
which satisfy the conditions: 

lim A = l, lim (B/r 2 ) = 1, lim C = 1. (97) 

r — >oo r — >oc r — >oc 

It allows to use all mathematical methods elaborated for general relativity. 
Conditions d96l l yield one equation: 

B 2 AC = r 4 . (98) 
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It allows to exclude the function A from the (|95l l. Then the equations Rn = and i?oo = are: 

-2BC' + 2rB'C + rBC" = 0, (99) 



-4BCB' + rCB' 2 - 2BC' + 2rBB'C + 2rBCB" + rB 2 C" = 0. (100) 

Because of equality ( [99b the sum of tree terms in ( llOOl i is equal to zero, and we obtain a differential 
equation for the function B: 

2rBB" + rB' 2 - ABB' = 0. (101) 

A general solution of this equation can be written as B = a(r 3 +/C 3 ) 2 / 3 where a and K, are some constants. 
The constant a = 1 which can be seen from condition d97l ). 

The function C can be found from differential equation (l99l in the form 



C" + 2 VB ' - C = 0, (102) 
rB 

where (rB' — B)/rB = (r 3 — /C 3 )/(r 4 + r/C 3 ). A general solution of this equation is C = b - Q/B 1 / 2 
where b and Q are constant. It follows from d97l ) that 6=1. 

Thus, in the spherical coordinate system the functions A, B and C are given by: 

A = 7HTWrY B = f ' c = 1 -T (103) 

where 

/ = (r 3 + /C 3 ) 1/3 . (104) 

V 

-a/3 



The nonzero components of the tensor Bl s are given by 



Rr- lA '- or _l 2rA-B\ \ 2rA-B' 

2 A ' tr ~ 2 C 6r ~ 2 rB ' * r ~ rB ' 

In non-relativistic limit the radial component of the equations of the motion of a test particle (Q]) takes 
the form x r = —c 2 Tq , where Tq Q = C /2A = r 4 C" / f 4 C. Therefore, to obtain the Newton gravity law 
it should be supposed that at large r the function f(r) coincides with r, and Q = r g = 2G M/c 2 is the 
classical Schwarzschild radius. 

At the given constant Q allowable solutions ( 1103b are obtained by changing the arbitrary constant JC. 
In particular, if we set JC = 0, then the line element (l94l l of the space-time in PFRs coincides with the 
Droste-Weyl solution of Einstein's equations [34] (it is commonly named Schwarzschild's solution) which 
has an event horizon at r = r g : 

fir 2 9 

ds 2 = ~ - r 2 [de 2 + sin 2 9 dip 2 } + (1 - rjr) dx° . (105) 

(1 - r g /r) 

If we set tC = Q, the line element coincides with the original Schwarzschild solution [35] 

f' 2 dr 2 9 

ds2 = - TT 77T - f^ d92 + sin2 9 + (1 - r g/f) > ( 106 ) 

(! - r a/f) 
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where / = ( r| + r 3 ) 1 / 3 . This solution has no event horizon and no singularity in the centre. Really, for 
motion of a particle in the plane 9 = it/2 laws of conservation of energy E and momentum J take place: 

&-&-L = E, d ± = J. (107) 

or otp oip 

It follows from this fact that the equations of the motion are of the form 

f 2 = (c 2 C/A)[l - (C/E 2 )(l + r 2 g J 2 /B)}, , (108) 



(f = c CJ/r g BE, (109) 

where (i, r, tp) are the spherical coordinates, f = dr/dt, ip = dip/dt , E = E/mc 2 , J = J/r g mc. 
The radial velocity of freely falling particle is given by the equation 



C) 



1/2 



(110) 



Fig.Q]shows the plot of the velocity as a function of the distance r/r g from the centre of a point attracting 
mass. It follows from this figure that this function is defined at all distances from the centre, and tends to 
zero at r — > 0. 




r / r 



Fig. 1 The velocity of a freely falling particle (in units of c) as the function of the distance r/r g from the centre. 



Line elements ( 1105b and ( 1 106b in V4 are not equivalent because they were obtained in the same co- 
ordinate system and correspond to the different values of the constant K, in solution ( 1103b . Both of the 
line elements refer to the spherical coordinate system, defined in £4, where t, r, tp and 6 are magnitudes 
measured by measurement instruments. Suppose, as it is usually believed, ( 11061 ) can be transformed to the 
form ( 1105b by an appropriate local coordinate transformation r — > r' and t — > t'. and [37]. But these 
new coordinates have not an operational meaning, and after an appropriate transformation of the physical 
3-space and time intervals, the same physical results will be obtained from d 1 05b as from ( 1106b . (Just as 
in classical electrodynamics in arbitrary coordinates). In particular, since the line element (1106b at K, = Q 
has no the event horizon and the singularity in the centre in spherical coordinates, it does not contain it in 
other coordinate systems. 



" Problems of correctness of such transformation have been considered in [34], and criticism with respect to the Droste-weyl 
solution in General Relativity — in [36] 
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It can be argued that the constant JC must be equal to Q = r g . Gauge-invariant tensor Bp (x) can 
be considered as the field strength tensor of gravitational field. There is a simple possibility to consider 
the object Bp (x) as a function of a tensor field ip a p(x) which may be interpreted as a potential i/j a p(x) 
described gravity in E4 as a field of spin two. Namely, one can set 

= V Q ^ 7 - (n + I)" 1 (<*|V CT Vv 7 + <5"V>/3 ff ) . 

The identity £>^ 7 = is satisfied as it is to be expected according to the definition of the tensor B^. Then, 
at the gauge condition V CT ^v 7 = Eqs. (l22l are 

= v^v^, v CT ^ 7 = o, (in) 

where □ is the covariant Dalamber operator in Minkowski space-time. 

It is natural to suppose that in the presence of matter these equations take the forrrQ 



Oip a p = x(T Q +t a p), 
V ct t/v 7 = 0, 



(112) 



where x = 8irG/c 4 , 



t a p = H~ X B U atl B» v = X^Wcry (H3) 

and T a p is the energy-momentum tensor of matter. 
Obviously, the equality 

V P {T aP + t a p) = (114) 

is valid. Therefore, the magnitude t a p can be interpreted as the energy-momentum tensor of the gravita- 
tional field. 

Let us find the energy of the gravitational field of a point mass M as an integral in Minkowski space- 
time: 



£gf = j hodV, (115) 
where dV is a volume element. In Newtonian theory this integral is divergent. In our case we have: 

t ao = 2x~ 1 B 1 m Bl 1 ^^-^- (116) 



and, therefore, in spherical coordinates, we obtain 

12 „4 



r 1 q 2 c 4 

£gf = J t ™ dV =8^G J > (117 > 



where [38] 
and 



I ^ = ^(1,1/3) (118) 



B(z,w)= r—dt (119) 

x ' 1 (1 + t) z+w 



It should be noted that, when we introduce ip a g in such a way, we cannot be sure a priori that the equations for ip a yield all 
physical solutions of the equations for B^. 
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is the B-function. Using the equality 

= WW 
1 (z + u>) 

where T is the r-function, we obtain J = 4n/IC, and, therefore, 

qqc 4 e 2 

^ = iC2G = iC Mc - (121) 

We arrive at the conclusion that at /C 7^ the energy of a point-mass is finite and, in particular, at 
K. = Q it is caused entirely by its gravitational field: 

£ gf =Mc 2 . (122) 

The spatial components of the vector P a = to a are equal to zero. 

Due to these facts we assume in the present paper that JC = Q = r g and consider the functions 

A = Pu^e/T) ' c = 1 ~ r ' /f - B = f2; (123) 

/ = (r 3 + , ,3)1/3 

in the Lagrangian d93l as a basis for analysis of physical effects. 

It must be noted that by using the above results we can find also the gravitational field inside a spherically- 
symmetric matter layer. In order to reach a coincidence of the equation of motion of test particles in nonrel- 
ativistic limit with the Newtonian one, the constant Q in ( 11031 ) in this case must be equal to zero. Therefore, 
the spherically-symmetric matter layer does not create gravitational field inside itself. This result will be 
used in next section. 

As an observer, located in a PFR does not feel influence of any forces, he should explain the fact of the 
deviation of paths of close particles of the PFR reference body as a manifestation of deviation of geodesic 
lines of these point masses in space-time with the line element ( 1106b . However, unlike the situation in 
General Relativity, he does not observe any problems at the Schwarzschild radius and close to r = 
because components of the curvature tensor have not any singularity. For example Rq 10 = r g f ~ 4 (r g — /). 
It tends to zero when r — ► 0. This fact is in accordance with the fact that in flat space-time the gravitational 
force tends to zero when the distance from the central point mass tends to zero, see next section. 



10 Cosmological test 

Physical consequences from the spherical-symmetric solution of equations ( f3Qb differ very little from those 
in general relativity at distances r ^> r g but they are different in principle at distances of the order of 
the Schwarzschild radius or less than that. Some principal physical consequences were considered in 
papers [40], [41], [42]. Some of them have chances to be checked up. The equations (f30b were successfully 
tested by the binary pulsar P5'i?1913 + 16 [39]. Here we consider another problem. 

Within last 8 years numerous data were obtained testifying that the most distant galaxies move away 
from us with acceleration [43]. This fact poses serious problems both for astrophysics and for fundamental 
physics [44]. In the present paper it is shown that available observed data are an inevitable consequence of 
properties of the gravitational force as deduced from the geodesic-invariant bimetric gravitation equations. 

According to (Q3, the gravitational force m f of a point mass M affecting a freely falling particle of 
mass m in Minkowski space-times is given by 

F = -m [c 2 C'/2A + (A'/2A - 2C' /2C)r 2 } , , (124) 

where the functions A and C are given by Eqs. dl23l ) and r is the distance from the centre. 
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For particles at rest (r = 0) 

GmM 



F 



1 



( r 3 +r 3)l/8 



(125) 



Fig. |2]shows the force F affecting particles at rest and the ones , free falling from infinity on a point mass, 
as the function of the distance r = r/r g from the centre. 




Fig. 2 The gravitational force (arbitrary units) affecting freely particles (curve 1) and particles at rest (curve 2) near 
an attractive point mass. 



It follows from this figure that the main peculiarity of the force acting on particles at rest in Minkowski 
space-time is that it tends to zero at r — > 0. The gravitational force acting on freely moving particles 
essentially differs from that acting on particles at rest. In this case the force changes its sign at ~ 2r g 
and becomes repulsive. Although we have so far never observed the motion of particles at distances of 
the order of r g , we can verify this result for very remote objects in the Universe at large cosmological 
redshifts, because it is well-known that the radius of the observed region of the Universe is of the order of 
its Schwarzschild radius. 

A magnitude which is related to observations of the expanding Universe is a velocity of distant star 
objects moving off from the observer. To find the radial velocity of such objects which is at the distance R 
from the observer, we can use spherically-symmetric solution ( 1123b of the vacuum equations d22l taking 
into account the effect of gravity of a spherically symmetric matter layer, as pointed out at the end of Sec. 
8. In this case M = (4/3)xpR 3 is the matter mass inside the sphere of the radius R, p is the observed 
matter density, and r g = (8/3)irc~ 2 GpR 3 is Schwarzschild's radius of the matter inside of the sphere. 

There is no necessity to demand the global spherical symmetry of matter outside of the sphere. Indeed, 
suppose, the matter density p at the moment is of the order of the observed density of the matter in the 
Universe (10~ 29 — 10~ 30 g cm -3 ). Fig. |3]shows the radial acceleration R of a particle in the expanding 
Universe on the surface of the sphere of the radius R. 

Two conclusions can be drawn from this figure. 

1. At some distance from the observer the relative acceleration changes its sign. If the R < 2 ■ 10 27 cm, 
radial acceleration of particles is negative. If R > 2 ■ 10 27 cm, it magnitute is positive. Hence, for suf- 
ficiently large distances the gravitational force affecting particles is repulsive and gives rise to a relative 
acceleration of particles. 

2. The gravitational force, affecting the particles, tends to zero when R tends to infinity. The same fact takes 
place as regards the force acting on particles in the case of a static matter . The reason of the fact is that at a 
sufficiently large distance R from the observer the Schwarzschild radius of the matter inside the sphere of 
radius R becomes of the same order as its radius. Approximately at R ~ 2 r g the gravitational force begins 
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Fig. 3 The acceleration of particles on the surface of a homogeneous sphere of radius R for three value of the 
parameter E. The matter density is equal to 10~ 29 ,g cm -3 . 



to decrease. The ratio R/r g tends to zero when R tends to infinity, and under these circumstances the 
gravitational force in the theory under consideration tends to zero. Therefore, the gravitational influence 
of very distant matter is negligible small. 

Thus, according to d 1 08b , the radial velocity of a star object at the distance R from the observer takes 
the form 



' R 2 



C_ 



(126) 



where C and / are functions of distance R of the object from an observer, v = dR/dt, E is the total energy 
of a particle divided by mc 2 . 

Proceeding from this equation we will find Hubble's diagram following mainly the method being used. in 
[45]. Let vq be a local frequency in the proper reference frame of a moving source at the distance R from 
the observer, v\ be this frequency in a local inertial frame, and v be the frequency as measured by the 
observer in the centre of the sphere. The redshift z = (v$ — v)jv is caused by both Doppler-effect and 
gravitational field. The Doppler-effect is a consequence of a difference between the local frequency of the 
source in inertial and comoving reference frame, and it is given by [22] 



= v Q [l- y/{l-v/c){l + v/c)]. 



(127) 



The gravitational redshift is caused by the matter inside the sphere of the radius R. It is a consequence 
of the energy conservation for photon. The energy integral (1107b for the motion of free particles together 
with the Lagrangian d93l l yield at v = and ip = the rest energy of a particle in gravitational field: 



E 



m c 



C. 



(128) 



Therefore, the difference in two local level energy E\ and E% of an atom in the field is AE = (E2 — 
E\)\[C, so that the local frequency v$ at the distance R from an observer are related with the observed 
frequency v by equality 



v = PiVC, , (129) 

where we take into account that for the observer location \[C = 1. It follows from ( 1129t and (1127b that 
the relationship between the frequency v as measured by the observer and the proper frequency vq of the 
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moving source in the gravitational field takes the form 



— = c 1 - V ^ C 

vq y 1 + v/c' 
which in the newtonian limit is of the form 



v 4 irGpR 2 N 



(130) 



(131) 



Equation ( 1130b yields the quantity z = v^jv — 1 as a function of R. By solving this equation numerically 
we obtain the dependence R = R(z) of the measured distance R as a function of the redshift. Therefore 
the distance modulus to a star object is given by 



fx = 5log w [R(z) (z + 1)] 



(132) 



where R(l + z) is a bolometric distance (in pc) to the object. 

If we demand that ( 1126l l has to give a correct radial velocity of distant star objects in the expansive 
Universe, it has to lead to the Hubble law at small distances R. At this condition the Schwarzschild radius 
r g = (8/3)irGpR 3 of the matter inside the sphere is very small compared with R. For this reason / r* r, 
and C = 1 — r g jr. Therefore, at E = 1, we obtain from d 1 26b that 



HR, 



where 



H = y/{8/S)wGp. 

If E 7^ 1 equation (1126b does not lead to the Hubble law since v does not tend to when R 



(133) 

(134) 
0. For this 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.! 



Fig. 4 The distance modulus fi vs. redshift z for the density p — 4.5 ■ 10 30 gem 3 . Small squares denote the 
observation data according to Riess et al. 

reason we set E = 1 and look for the value of the density at which a good accordance with observation 
data can be obtained. Fig. d4]i show the Hubble diagram obtained by (1132b compared with the Riess et al 
data [43] . It follows from this figure that the model under consideration does not contradict observation 
data. With the value of the density p = 4.5 • 10" 



~ 30 <7 cm 3 we obtain from ( 1134b that 



H 



1.59 ■ lO-^c" 1 = 49 (km/s)/Mpc. (135) 

Fig. |5]shows the dependence of the radial velocity v on the redshift. It follows from this figure that at 
z > 1 the Universe expands with an acceleration. At R — > oo the velocity and acceleration tend to zero. 
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Fig. 5 The radial velocity vs. redshift z for the density p — 4.5 ■ 10 30 g cm 



11 Conclusion 

The necessity of a discussion of the physical description of gravity through geodesic-invariant equations 
is quite obvious. However, the satisfactory implementation of such a programme is a difficult task. The 
equations of gravitation considered in the given paper do not contradict observation data, and their some 
predictions seem rather attractive. However it is still unclear whether they are only possible equations. 
Besides, it appears that the correct equations of this type should be the bimetric ones. However physical 
interpretation of such bimetricity involves deep problems of the space-time theory which presently have 
no a satisfactory understanding and a complete solution. 
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